Interference of a Bose-Einstein condensate in a hard-wall trap: 

Formation of vorticity 
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We theoretically study the coherent expansion of a Bose- 
Einstein condensate in the presence of a confining impene- 
trable hard-wall potential. The nonlinear dynamics of the 
macroscopically coherent matter field results in rich and com- 
plex spatio-temporal interference patterns demonstrating the 
formation of vorticity and solitonlike structures, and the frag- 
mentation of the condensate into coherently coupled pieces. 
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A landmark experiment |l) demonstrated in a strik- 
ing way the interference of two freely expanding Bose- 
Einstein condensates (BECs). In this paper we theo- 
retically study the evolution of a BEC in a coherently 
reflecting hard-wall trap. The present situation is closely 
related to the recent experiments on an expanding BEC 
in an optically-induced 'box' by Ertmer and coworkers 
1^ . Due to the macroscopic quantum coherence the BEC 
exhibits rich and complex self-interference patterns. We 
identify the formation of vorticity and solitonlike struc- 
tures, and the dramatic fragmentation of an initially uni- 
form parabolic BEC into coherently coupled pieces. 

Atomic BECs exhibit a macroscopic quantum coher- 
ence in an analogy to the optical coherence of lasers. 
In the conventional reasoning the coherence of a BEC 
is introduced in the spontaneous symmetry breaking. 
Nevertheless, even two BECs with no phase informa- 
tion, could show relative phase correlations as a result of 
the back-action of quantum measurement Moreover, 
the density-dependent self-interaction of a BEC demon- 
strates the analogy between nonlinear laser optics and 
nonlinear atom optics with BECs. BECs are pre- 
dicted to exhibit dramatic coherence properties: The 
macroscopic coherent quantum tunneling |6| and the for- 
mation of fundamental structures, e.g., vortices [0-^ and 
solitons UJl^. Some basic properties of grey solitons 
have been recently addressed for harmonically trapped 
ID BECs in Ref. [|o). Also optical solitons have been 
actively studied in the ID homogeneous space |Q . 

In this paper we study the dynamics of a BEC con- 
fined to a hard- wall trap with potential ^(r). Such walls 
can be realized, e.g., with a blue-detuned far-off-resonant 
light sheet. Throughout the paper we focus on repulsive 
interactions. When the BEC is released from a magneto- 
optical trap (MOT) inside the confining potential by sud- 
denly turning off the MOT, the repulsive mean-field en- 
ergy of the condensate transforms into kinetic energy and 
the BEC rapidly expands towards the walls. The reflec- 



tions of the matter wave from the binding potential re- 
sult in rich and complex spatio-temporal and interference 
patterns referred to in ID as quantum carpets p^ . They 
have been recently proposed as a thermometer for mea- 
suring the temperature of the BEC |13|. In this paper 
we show that the nonlinear dynamics of a BEC displays 
dramatic local variations of the condensate phase includ- 
ing the formation of vorticity and solitonlike structures. 
The solitary waves could possibly be used as an experi- 
mental realization of the macroscopic coherent tunneling 
analogous to the Josephson effect. 

The dynamics of a BEC follows from the Gross- 
Pitaevskii equation (GPE) 
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Here M and k = iTrh^aN/M denote the atomic mass 
and the coefficient of the nonlinearity, respectively, with 
the scattering length a and the number N of BEC atoms. 
Our initial distribution ip{r;t = 0) is the stationary solu- 
tion of the GPE with the potential V{r) replaced by the 
potential of the MOT. 

We integrate GPE in one and two spatial dimensions. 
The projections from 3D into ID or 2D require that the 
mean field ip does not vary significantly as a function 
of time in the corresponding orthogonal directions. This 
condition can be satisfied, e.g., in the presence of a strong 
spatial confinement to these dimensions. Then we can ap- 
proximate the position dependence in these directions by 
constants A and £ resulting in V'(r) — iJi{x)/A^^^ in ID 
and ■0(1") ~ "02(2^, y)/^^^^ in 2D. This yields the strengths 
Ki — k/A and K2 — k/I of the nonlinearity in GPE for 
the mean fields in ID and 2D ipi and ip2- We emphasize 
that especially the 2D calculations may already contain 
the essential features of the full 3D coupling between the 
different spatial dimensions by the nonlinearity. 

One- dimensional case - The linear Schrodinger equa- 
tion for the ID box of length L exhibits ||l^ regu- 
lar spatio-temporal patterns. These patterns consist of 
straight lines, so called traces, of different steepness corre- 
sponding to harmonics of a fundamental velocity, vq. The 
traces arise from the interferences between degenerate 
eigenmodes of the system. The eigenmodes of frequency 
ujn = n?LOi = n^hn"^ /{2ML'^) are a superposition of right 
and left propagating plane waves with wavenumbers fc„ = 
WK / L = nki. The probability density consists of the in- 
terferences between different eigenmodes. Therefore the 
lines of constant phase ±{km^kn)x + {ijJm —^n)t = const, 
correspond to straight lines in space-time with velocities 
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Vmn = ±(wm " w„)/(fcm ± = (m ± n)vo. Here we in- 
troduced the fundamental trace velocity vq = hn/ {2mL). 
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FIG. 1. [Part (a) in a separate file] Quantum carpet of a 
BEC woven by propagating grey solitonlike structures. Time 
evolution (a) of a BEC density profile |'!/;i(a;; in a ID 
hard- wall trap. The horizontal and vertical axes denote the 
spatial coordinate and time, respectively. The initially con- 
fined BEC expands towards, and reflects from, the hard-wall 
boundaries. The dark canals correspond to low density and 
represent the evolution of grey solitonlike structures. The 
phase slip of a single fringe (b) at t/T ~ 3.2 x 10~^ , and of four 
colliding fringes (c) at t/T ~ 0.021. The arrows in (a) denote 
the times of the phase profiles (b) and (c) . Here time is scaled 
in units of the revival time T = 2-k / ui = AML^ /{-kK) of the 
linear Schrodinger equation and the initial state is R/ L ~ 0.3 
and = 19000ft/T. 

We now turn to the quantum carpet of Fig. ^ rep- 
resenting a ID BEC trapped between two impenetrable 
steep Gaussian potentials approximating infinitely high 
walls dX X = ±L/2. At time t = the BEC is released 
from a harmonic trap of frequency fi. In the limit of 
strong confinement the initial state is well approximated 
by the Thomas- Fermi solution il}i{x\t = 0) = d{Ri ~ 
\x\)[i{Rl - x')/{2Rl)Y/^. Here i?i = [3Ki/(A/r!2)]i/3 
describes the ID radius of the BEC wave function. 

After the turn-off of the MOT the kinetic energy term 
in GPE becomes dominant and the matter wave expands 
towards, and eventually reflects from, the box bound- 
aries. Due to the macroscopic quantum coherence of a 
BEC different spatial regions of the matter field generate 
a complex self-interference pattern that exhibits canals 
analogously to the quantum carpet structures of the lin- 
ear Schrodinger equation ||l2|. From Fig. ^we note that 
only destructive interference fringes, canals, appear in 
the nonlinear carpet. Constructive interference fringes, 
ridges, do not emerge. 

For the single particle Schrodinger equation the result- 
ing quantum carpet demonstrates the fundamental wave 
nature of the particle. Therefore, it is perhaps surprising 
that in the case of GPE, which represents the coherent 
matter field of the many-particle system, the intermode 
traces acquire properties that demonstrate dramatic par- 
ticle nature. In particular, in the case of repulsive inter- 
particle interactions the traces correspond to the evolu- 
tion of solitonlike structures with an associated phase 
kink in Fig. |l|b,c. 

Grey solitons correspond to the propagation of 'density 
holes' in the matter field and the depth of the hole charac- 
terizes the greyness of the soliton. In homogeneous space 



the soliton is called dark, when it exhibits [g_l|jlC|] a van- 
ishing density at the center of the dip and a sharp phase 
slip of TT. For grey solitons the value of the phase slip 
Lp is reduced, thereby giving the soliton a nonvanishing 
speed according to v — ccos{ip/2), where the maximum 
velocity c = {inh^ap/M^y/^ is the Bogoliubov speed of 
sound at a constant atom density p. As a result of the for- 
mation of the soliton-like structures the initially uniform 
parabolic BEC is dramatically fragmented into spatially 
separated pieces and we can interprete these density holes 
as the boundaries of the fragmented BEC. Analogously 
to the Josephson junction the phase of the BEC is ap- 
proximately constant outside the narrow region of the 
solitary wave. The coherent tunneling of atoms across 
the boundary depends on the relative phase between the 
two contiguous pieces and results in the motion of the 
'solitary wave junction'. For a phase slip of tt the veloc- 
ity of the soliton vanishes analogously to the vanishing 
Josephson oscillations of the number of atoms. 

The interatomic interactions shift the eigenenergies 
compared to the linear case. As a result the fringe ve- 
locities in Fig. are a few tens of percents larger than 
the linear trace velocities Auvq for the canals with a sym- 
metric initial state p^ . The nonlinear fringe velocities 
are still determined by the degenerate intermode traces 
and are approximately integer multiples of the minimum 
speed. The largest velocity in Fig. ^ is approximately 
given by the Bogoliubov speed of sound in terms of the 
average density of atoms in the box. 

The fringe evolution shown in Fig. |l| displays a remark- 
able inherent particle-like robustness and solitonlike be- 
haviour: The fringes survive complex dynamics. Their 
paths exhibit dramatic avoided crossings demonstrating 
repulsive interactions and therefore the colliding wave 
packet holes become degenerate in momentum space . 

Two-dimensional case - We now turn to the evolu- 
tion of a BEC in a circular box. Again we approximate 
the infinitely high walls by a steep Gaussian potential 
aligned along a circle of radius r = L/2. The initial wave 
function is the Thomas-Fermi solution of the symmetric 
MOT ^P2{x,y;t = 0) = [2{Rl - x"^ - y'^)/{nR^)Y/'^ , and 
i?2 = [4K2/(7rMf]2)]i/4 denotes the 2D 'classical' radius. 
Hence, the initial wave function is located at the center 
of the circular box and we can show that in this case the 
state remains symmetric also at later times. 

Figure ^ shows the 2D density |-!/'2(a;, y; i)P and the 
phase profile \(p\ of a BEC obtained from GPE at a later 
time. We note the formation of a regular interference 
patterns similar to solitary waves. The fragmentation 
of the BEC is even more dramatic than in the ID case: 
The fringes exhibit a vanishing density at the center of 
the dip. The BEC forms coherently coupled loops and 
the resulting structures are similar to optical ring soli- 
tons 

When we now slightly displace the initial state of the 
BEC from the center of the circular box, the rotational 
symmetry is broken. In Figs. ^ and ^ we show the result- 
ing 2D density (left column) and the phase profiles (right 
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column) at three characteristic times. The reflections of 
the BEC from the hard-wall potential create solitonlike 
structures. At later times the stripes bend and eventually 
break up forming dark spots that correspond to vortices 
with associated phase windings around closed paths wj. 




FIG. 2. Formation of loops of solitonlike structures in a 
BEC expanding in a 2D circular box. The symmetric con- 
densate has started from the center of the circular box, has 
expanded and reflected from the boundary. We show the 
density (left) and phase (right) profiles of the BEC at time 
t = 4.9 X 10~^ M /h. In the phase profile we have chosen a 
continuous phase \(p\, where — tt < < tv. The initial state is 
R2/L ~ 0.28 and K2 = 1500h^ /{ML'^). 

We note that we can use the present situation of an 
expanding BEC in a circular box to create vortices in 
some particular spatial location by simply introducing a 
static potential dip and letting the expanding BEC flow 
across it. This is a simplified version of the suggestion by 
Ref. Q that a moving potential barrier through a BEC 
can create vorticity in the vicinity of the potential. 

As a final example and to demonstrate the effect of 
the symmetry of the hard- wall trap we consider the evo- 
lution of a BEC in the 2D square box. We generate the 
boundary by steep Gaussians approximating infinitely 
high walls a,t x = zLL/2 and y = ±L/2. 

The square has the symmetry of rotations of 7r/2. 
Hence, for a symmetric nonrotating initial state the min- 
imum number of vortices conserving the total angular 
momentum is eight. In Fig. |^ we show the density pro- 
file at two characteristic times. The reflections of the 
BEC from the boundary generate dark solitonlike struc- 
tures with an amazingly regular square shape that start 
bending, break up, and form vorticity. 

We now address the experimental feasibility of the 
proposed self-interference measurements. A higher-order 
Laguerre-Gaussian beam can generate a hollow optical 
beam Jist and thus a cylindrical atom-optical hard-wall 
potential around the magnetically-trapped BEC. More- 
over, BECs with purely 2D confinement may also be in- 
vestigated in certain magnetic trap configurations 
In recent experiments Bongs et al. studied the coher- 
ent reflections of a BEC by atom-optical mirrors and the 
evolution of a BEC in a atom-optical waveguide closely 
related to the present theoretical study. 



FIG. 3. [In a separate file] Formation of vorticity in a 
BEC expanding in a 2D circular box with displaced initial 
state. We show the density profiles (left column) and the cor- 
responding phase profiles (right column) of the BEC at three 
characteristic times t = 2.5x 10"^, 4.9 x 10~^, and 6.1 x 10"^ 
in units of ML^ /h. The reflections from the boundary gener- 
ate sofitonlike structures that bend, break up, and form vor- 
ticity. The velocity of the BEC is proportional to the gradient 
of the phase. 

The density profiles of the BECs could be directly mea- 
sured via absorption imaging, if the necessary spatial res- 
olution could be obtained, e.g., via ballistic expansion of 
the atomic cloud. Vortices may also be detected by inter- 
fering a BEC with and without vorticity. Then the phase 
slip in the interference fringes would be the signature of 
the vorticity [^. The phase slip between two pieces of 
the BEC has a dramatic effect on the dynamical struc- 
ture factor of the two-component system w hich may 
be observed, e.g., via the Bragg scattering ([iSj. 

In conclusion, we studied the generation of vortic- 
ity and solitonlike structures of a BEC in a hard-wall 
trap. The nonlinear evolution of GPE dramatically di- 
vides the initially uniform parabolic BEC into coherently 
coupled pieces. We showed that the density profile of 
the BEC can be a direct manifestation of the macro- 
scopic quantum coherence. Obviously it could also be a 
sensitive measure for the decoherence rate of the BEC 
JlOt . Unlike the typical coherence measurement that de- 
tects the relative macroscopic phase between two well- 
distinguishable BECs |p],p^, the present set-up probes 
the self-interference of an initially uniform matter field. 

We acknowledge financial support from DFG and from 
the EC through the TMR Network ERBFMRXCT96- 
0066. 




FIG. 4. The magnification of the last density and phase 
profile of Fig. ^ displaying vorticity. For example, the dark 
spot in the center represents a vortex with a unit circular 
quantization of 2n. 
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FIG. 5. Expanding BEG in a 2D square box with a 
symmetric initial state at the center of the square. At 
t = 2.9 X W'-^ML^/h (left) the density profile of the BEG 
exhibits a rectangular pattern of solitonlike structures. At a 
later time t = 3.9 x 10~^ML'^/h (right) the pattern is dis- 
torted and we see the beginning of the formation of vorticity 
symmetrically around the diagonals. The initial radius of the 
BEG is R2/L ~ 0.28 and K2 = 5000h'^ /{ML'^). 
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